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Abstract 



We study a fermionic infinited-ranged Ising spin glass with a 
real space BCS interaction in the presence of an applied trans- 
verse field. The problem is formulated in the integral functional 
formalism where the SU(2) spins are given in terms of bilinear 
combinations of Grassmann fields. The problem is solved within 
static approximation and the replica symmetry ansatz combined 
with previous approaches used to study the critical behavior of 
the quantum Ising spin glass in a transverse field and the spin 
glass Heisenberg model with BCS pairing. Our results show that 
the transverse field has strong effect in the phase boundary of 
the spin glass phase and the PAIR phase in which there is a long 
range order corresponding to formations of pairs. The location 
of the tricritical point in the PAIR phase transition line is also 
affected. 



PACS numbers: 05.50. +q, 6460.Cn 



1 Introduction 



Theoretical studies in recent years have been investigating the interplay be- 
tween superconductivity (SC) and spin glass (SG) [H 12 El HI E| which has 
been found in several strongly correlated electron systems, such as heavy 
fermions [Oj and cuprate superconductors [7 j . However, the experimental ev- 
idences for these correlated systems have showed a quite complex situation. 
For instance, the heavy fermion superconductor Ui- x M x Pd,2Al^ (M = La, Y, 
Th) [6 j shows a sequence of magnetic and Non-Fermi Liquid (NFL) ground 
states. In particular, when the content of La is increased, a antiferromag- 
netic phase is replaced by a SG which is suppressed to T = OK at x — 0.75. 
After the Quantum Critical Point (QCP), there is NFL region and at x — 1 
appears superconductivity. 

Some approaches have proposed that disorder itself can be the source of 
deviation of the Fermi Liquid behavior. For heavy fermions, for instance, 
the so called Kondo disordered model (KDM) jS] describes a distribution of 
Kondo temperature Tk given origin to the NFL behaviour. Castro-Neto et 
al. |H] relates the NFL effects to the presence of an inhomogeneous Griffith's 
phase. In an earlier work, it has been shown a NFL behavior near to the 
QCP. in a transition between a metallic-paramagnet and a metallic-spin glass 
[TP] . But, there is also the suggestion that the presence of QCP itself can be 
a source of a new class of excitations which spread its effects even at finite 
temperatures leading to a breakdown of the Fermi liquid theory [TT] . 

However, little consideration has been given to describe how the phase 
boundaries between superconductivity and spin glass are modified when at 
low temperature, besides the presence of thermal fluctuations, the quantum 
fluctuations start to become important. Some of the models previously men- 
tioned [21 m h ave addressed the phase transition problem between SG and 
superconductivity. Nevertheless, in those references there is no mechanism 
able to tune quantum fluctuations. 

In particular, the model in Ref. [2] has been derived from a model in- 
troduced by Nass et al |T2] to deal with conventional superconductors doped 
with magnetic impurities. This model is an s-d exchange interaction be- 
tween the magnetic impurities together with a conventional BCS interaction 
between the conducting electrons. When the conduction electrons are in- 
tegrated by a perturbation expansion to second order in J sd (the exchange 
interaction), the resulting effective model consists of the RKKY interaction 
with a pairing interaction between localized fermions. If the coupling 
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between the localized magnetic moments is assumed to be random gaussian 
distributed, one has a model to study the phase boundaries between SG and 
a phase where there is spin pairing J2J . 

The functional integral machinery, where the SU(2) spins have been rep- 
resented by bilinear combination of the Grassmann variables, has been the 
suitable method used [2] to study the mentioned effective model in its Ising 
version. Combined with replica trick and the static approximation, it has 
allowed to find the Grand- Canonical potential in terms of the spin glass and 
the PAIR order parameters. The last phase corresponds to a long range order 
where there is pair formation. In the half-filling situation, the results have 
showed a phase diagram temperature T/J versus g/J (g is the strength of 
the pairing interaction and J is the variance of J^) where one can find a spin 
glass phase for low temperature and small g. If g is increased, one gets a 
phase transition at g — gi(T) where there is a PAIR phase. The nature of 
the transition line is complex presenting a tricritical point (T ic , g tc ). 

Further investigation studied the pairing-spin glass competition replacing 
the fermionic Ising in the cited effective model by the Heisenberg model with 
an applied magnetic field H z within the same framework of References j2] 
and [13J. The results have showed that the transition to the PAIR phase de- 
pends on the replica diagonal spin glass order parameter (which is associated 
with the susceptibility) even at higher temperatures than the freezing tem- 
perature Tf. The region in temperature where the calculated line transition 
between the normal-paramagnetic (NORMAL) and the PAIR phases is first 
order becomes larger when H z is increased. Nevertheless, the line transition 
between the SG and the PAIR phase have not been accessible in this work. 
Other interesting point about the infinite ranged quantum Heisenberg model 
in the presence of a field is that Tf is depressed when the field increases, but 
never reaches a QCP [3]. 

Corrections in the weak hopping limit to the Ising SG fermionic model 
in the presence of a local pairing interaction have been studied elsewhere [3]. 
The results show that those corrections essentially preserves the shape of the 
phase diagram obtained with no hopping. Therefore, one can consider the 
the range of validity of this theory as covering the transition between poor 
conductors and superconductors. 

Recently, the Ising SG alone has been investigated in the transverse field 
T. The functional integral approximation have been used ^3] to deal with 
the non-commutativity of the spins operators which have been represented 
by bilinear combinations of Grassmann fields. There are two versions of SG 
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fermionic problem. The first one, the operator Sf has four eigenvalues (two of 
them are non- magnetic). In the second model, the two vanishing eigenvalues 
are suppressed by a restraint. In both models, the freezing temperature de- 
creases with increasing V until to reach a QCP at a critical value T c This sort 
of approach is a natural tool to study phase transitions in condensed mat- 
ter problems where fermions experiment couplings such as superconductivity 
and Kondo effect |To"j . 

Therefore, our aim in this work have been to investigate how the phase 
boundary between a PAIR phase (where there is pair formation) and the four- 
state fermionic spin glass (SG) is modified if there is spin flipping induced by 
the presence of the transverse field T. Therefore, the transverse field plays 
the role of the spin flipping part of Heisenberg model [TH] allowing to access 
the QCP. In order to solve the functional integral over the Grassmann fields 
contained in the partition function, the formalism of Nambu matrices and 
spinors has been used. We also have used the replica symmetric "ansatz", 
and therefore we have calculated the Almeida-Thouless line to obtain 
the validity limit of this procedure. Finally, we find the Grand-Canonical 
potential and the saddle point equations for the order parameters in the 
half-filling limit. 

One important approximation in the present work is the neglecting of time 
fluctuations (the static approximation) [T£]. For infinite ranged quantum 
Ising spin glass, a simulational approach jTU] has shown that for temperatures 
close to the freezing temperature Tf (when T = 0) the static approximation 
can be considered reliable. On the other hand, it is quite clear that the 
static approximation is unable to capture the fundamental low temperature 
dynamical behavior of the correlation functions [2D] ■ Nevertheless, it has been 
shown that for M-component quantum rotor model, in the limit M = oo, 
the critical line is given by zero-frequency mode [21] . This critical behaviour 
coincides with the Ising SG in transverse field [20]- Therefore, that is the 
ultimate justification for the use of the static approximation to find the phase 
boundary between SG and the PAIR phase for increasing transverse field 
which is the main purpose of the present work. 

This paper is organized as follow. In section II, we introduce the model 
and perform calculations using the replica trick and the static approximation 
in order to find the Grand canonical potential and the saddle point order 
parameters equations. The behavior of the tricritical point in the transition 
line to the PAIR phase is obtained as a function of both T and g. In section 
III, phase diagrams are build up with solutions from the set of the order 
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parameter equations in both situations, T/J versus g/J (for two values of 
T) and T/J versus T/J (for several values of g). It is also suggested a 
relationship between V and g which allow to see more clearly the role of 
quantum fluctuations in the interplay between SG and the PAIR phase. In 
the last section, we present our conclusions and final remarks. 



2 General Formulation 

The model considered in this work was obtained by tracing out the conduct- 
ing electrons degrees of freedom of a superconductor alloy j2], resulting in 
an effective BCS pairing interaction among fermions and a random Gaussian 
interaction coupling the localized spins. In the resulting effective model we 
apply the transverse field term. Therefore, the Hamiltonian is 



jijS-Sj + ^4T c !|Cjj.Cjt 

-£2r£? + J> £ n js (1) 

3 3 9=1,1 

where the sum is over the N sites of a lattice. The coupling is an inde- 
pendent random variable with Gaussian probability distribution given by 



The spin operators in Eq. (JTJ) are defined (see references [312]) as: 



1 1 

S j = - fyi] ; S j = ^[sVii + c \i c 3l\ ( 3 ) 



where the c\„ ( fermions creation (destruction) operators, w ith a =] 

or I indicating the spin projections, fij a = Cj a Cj c is the number operator and 
/i is the chemical potential. 

The Grand Canonical partition function is formulated in the functional 
integral formalism for fermions using the anticommuting Grassmann variable 
0j s (r) and 0j S ( r )( T is the complex time). Therefore one has 

Z = J D{<f)*(f)) exp [A + A SG + A r + A BCS ] (4) 
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where the actions A , A SG} A r , and A BC s are the free part, the spin glass 
part, the transverse field part and the pairing part, respectively. The three 
first ones assume (after time Fourier transformation) the following forms: 



4 = EEi t M( iw + ^)4 

j u> 

A SG = J2i3j ij sms^-n), 

ij 



(5) 
(6) 
(7) 



where S?(Q) = |X^0 .(u; + Q)g_ 3 4> . (u), with Matsubara's frequencies uj 



(2m + 1)tt and Q = 2m7i (m = 0, ±1, ■ ■ •). In the Eqs. ©-© we have used 
the Spinors 



0>) 



and the Pauli matrices 



1 

1 



Ol 



-i 

1 



1 
-1 



The pairing action is given by 



A BCS = T,T,pmpM 

ij Q 



(8) 



(9) 



(10) 



With Pj(n) = Eo; 0j|(-^)0jT( fi + <*>)■ 

In this paper we discuss the phase transition problem within the static 
approximation, therefore only the term with Q = is kept in the sum over 
the Matsubara's frequencies in Eqs. (JUJ) and (jlO)) . For this reason, we can 
define Sj = 5*1(0) and write the following expression for the pairing action 



.4 



St 

~BCS 



£g_ 

AN 



p=l,2 



3,u 



(11) 



where the Nambu matrices have been introduced in the previous equation, 



0>) 



(12) 
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The configurational averaged Grand Canonical potential per site can be 
found by using the replica formalism 



n 



1 Z(n)-1 
lim 



N Nf3 n^o n 



(13) 



The configurational averaged replicated partition function Z(n) = (Z n )j t , 
becomes, after averaging over J^-, 

Z(n) = J « exp j^^g + A?]+ 

a p=l, 2 \j> / 



?2j 2 



AT 



(14) 



a,/3 \ 3 



In the previous equation we introduce the replica index a = 1, 2, - • • , n. 
The linearization of Eq. is obtained by using the Hubbard-Stratonovich 
transformation 

/oo />oo ^ 

Y[driRadr) Ia / exp {-iV 



a/3 



X 



/3 2 J 2 



lip - lnA Q (g a/3 ,?7 a ))}, 



(15) 



tt/3 



where was introduced replica dependent auxiliary fields \q a \ and q a p. In the 
Eq. 7?., = - A/ = (0gN/n) n (0 2 J 2 N/2'n;) n , where 



A a (g a/3 ,?7 C 



a=l 



exp 



+4/? 2 j 2 x: + e^»zl 

a/3 w 



(16) 



with the matrix r? defined as: 

_Lrv 



/^^r/o 
/3gri* a 



(17) 
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The order parameter \r] a \ introduced in the Eq. (|15|) corresponds to a long 
range order where there is double occupation of the sites |2], and q a p is the 
spin glass order parameter. 

In the present work we restrict the discussion to the replica symmetric 
ansatz, that considers 



q a p = q; q c 



q + x 



(18) 



where q is the spin glass order parameter and % — § (x * s the static sus- 
ceptibility [22j). One can sum over the replica indices, which produces 
new quadratic terms that are linearized by introducing new auxiliary fields. 
Therefore, the functional integral becomes 



Dz 



(19) 



with Dz = dz^^, = dt^S^- , and 



2tt 



2tt 



I{(,z,h) = jl)[0*0]exp 



OJ 



(20) 



where the matrix G x (uj) is defined by G_ x (uj n ) = iu n + + PTa± + ha 3 
and the field h = fiJ^WZ + (3J^2qz. 

In order to solve the integral in Eq. (J2"U|) which combines the elements of 
spinors and Nambu matrices, we can use a similar procedure already done in 
Ref. ||| which mixes the elements of the spinors and the Nambu matrices to 
write Eq. fl2U as: 



I(S,z,h) = j D[<f0]exp 



(21) 



where 



(22) 
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and 



( iuj + C (3T Pgr} \ 

PT + C -Pgr} 

Pgr}* ico-C + ~PT 




(23) 



\ u -pgr}* -pT J 

with C ± = Pfi ± h. 

In the Eq. (J2TJ), the differential £>[</>*</>] stands for FL Ua=n ^C(^) 
d</>* (— cj) dcpa-i^uj) d(f) a (—L)). The functional integral over the Grassmann fields 
and the sum over the Matsubara frequencies can be readily performed with 
the result 013]: 



J(f , *, h) = cosh J^/i) 2 + (Pg\r}\) 2 + cosh V6 



(24) 



where 6 = h 2 + (PT) 2 . 

The results obtained in Eq. ()24)1 can be used in the Eq. that allow us 
to rewrite the Eq. (|15|) . Therefore, the saddle point method (see Eq. (|13|) ) 
give us the Grand Canonical potential as 



where 



~- P9V 2 + ^X(X + 2q)- Dz In I a (z) 

L J — oo 

/oo 
Dtl cosh v9 
-oo 



(25) 



(26) 



with the chemical potential fixed to ensure that we are in the half-filling 
situation. From now, the parameter rj is used instead of \rj\. The saddle 
point equations for order parameters that follow from equation (J25|) are: 



Dz 



<-™ swh(Pgri) 
2J-oo Z Uz) 



fZoDth sinhv/6 



Uz) 



X 



-OO 



Uz) 



(27) 

(28) 
(29) 
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where 

'^=L D( (^-ij — + )■ (30) 

The solution with replica symmetric is unstable at low temperature. 
Therefore, it is necessary to find the region in temperature where the Almeida- 
Thouless eigenvalue [T?] Xat becomes negative, which is given by: 



r2' 



2 



A w = 1 - 2( >J) 2 r Dz ( I ^f§ ) h ) ■ (31) 



In the Eq. (JHU I c (z) is 



W) = r (32) 



The Landau expansion of the Grand Canonical potential in powers of the 
two parameters r] and q allow us to locate the second order line transition in 
the problem. Therefore, from Eqs. (j25|1 - (j2f)j) : 

pn = j2fMx,P)<i j (33) 

3=0 

where x(g, rj, (3) is the solution of the saddle point equation, that was rewrite 
in powers of q as: 

X = Xo + XiQ + X2q 2 ■ (34) 

Introducing the Eq. into Eq. (J3~3*j) . and expanding the coefficients fj 
in powers of q and 77, we obtain the following result: 

^ = ^fxl ~ HKo) + A 2 q 2 + A 3 q* 

+B 2V 2 + B 4 r) 4 + B 6V 6 (35) 



with 



A 2 = -^f + ^x 2 , A 3 = - 8 -^xl (36) 
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Figure 1: Phase diagrams as a function of T / J and pairing coupling g/J for 
two values of V. Solid lines indicate second-order transition while dotted line 
indicate a first-order transition. The tricritical point for T = is shown in 
detail in the diagram. 



Bo 



(3g 



Ba 



P 4 g 



A\K 2 



gv 

6!K n 3 



/°° D^(coshv / A-2), 

J — oo 

3—13 Z" 00 Df coshVA 

«/ — OO 

+ (/°° £>f coshv'A 



(37) 
(38) 



where -Kn 



(39) 

1 + H^coshVA , ^ = MWA , an d A = 

2p 2 J 2 M 2 + P 2 r 2 . 

The tricritical point is given when both coefficients B2 and B4 change the 
sign. In this condition we have from Eqs. (}3*Tj) and (}3*Hj) that: 



/oo 1 
D£ cosh y/A t = 2, 
-00 



(40) 
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/3 t gt = 6 



(41) 



where A t is defined similar to A above with: 

1 r°° sinhi/A^ , , 

where the sub-index t stands for the tricritical values of the temperature T, 
g and V. We have solved numerically Eqs. (f4~T|) and (|42|) for several values of 
T t . The result are shown in the phase diagrams (see Figs. (JTJ), (J2J) and (3)). 



3 Results 

The numerical solutions of Eqs. (}2T)) -(f2T H) allow one to construct two sorts 
of phase diagrams when T and g are independent parameters. The first one 
is T/J (T is the temperature) versus g/J (g is the strength of the pairing 
interaction) where the transverse field T/J is kept constant. The second one 
is T/J versus T/J with g being constant. 

In the Fig. 1, one can see the results obtained in a phase diagram T/J 
versus g/J for T/J = and T/J = 1 (for numerical purposes J = 1). In 
the first case, we have obtained the same phase diagram already found in 
Ref. |2j with three distinct regions. The normal-paramagnetic (NORMAL) 
region at high T and small g (where q = and r] = 0). For g > gi(T) 
(see section 1) one enters in the PAIR phase (where q = and rj ^ 0). 
Finally, for low T and small g, one has the phase transition to the spin glass 
phase (SG) at T = Tf (Tf is the freezing temperature). When T is turned 
on, the freezing temperature decreases and the line transition g = g±(T) 
is displaced showing a dependence with T. Therefore, it is necessary to 
increase simultaneously the parameter g to find again solutions of the order 
parameters (see Eqs. (J2Tj) - (f2§j) ). which corresponds to the PAIR phase. The 
position of the tricritical point (T tc , g tc ) also moves when T is increased, and 
the multiple solutions region corresponding to the first order line transition 
is enhanced. We also have obtained the behavior of the Almeida-Thouless 
(AT) eigenvalue Xat showing that for both values of T, the replica symmetric 
SG solution is unstable. 

In Fig. 2, the phase diagram is plotted T / J against T / J for several values 
of g. For g = (see Fig 2. a), the corresponding phase diagram reproduces 
basically the results found in Ref. jT5j. These results show that the freezing 
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g = O.OJ (a) 
NORMAL PHASE 


g = 6.5J (b) 

- — NORMAL PHASE 
PAIR 

phase h^--. 




— SG PHASE 


- _ g = 8.0J (c) 

\ NORMAL PHASE 
PAIR PHASE ) 

>-M 1' 

: " SG PHASE ^\ 


^\^=9.0J (d) 

'-.NORMAL PHASE" 

PAIR PHASE 
- SG PHASE 



0.5 1 1.5 2 2.5 r 0.5 1 1.5 2 2.5 r 

c c 

r/j r/j 

Figure 2: Phase diagrams as a function of T/J and Y/J for several fixed 
values of g/J: (a) g = 0, (b) g = 6. 5 J, (c) g = 8 J, and (d) g = 9 J. It is used 
the same convention as Fig. for the transition lines. The hatched region 
in panel (b) delineates a multiple solution region. 

temperature Tf decreases (when Y increases) towards to a QCP with r c = 
2^/2. The entire SG region in the phase diagram is unstable (see the A-T 
line in Fig. (2. a)). If g is turned on, which energetically favors the double 
occupation, the PAIR phase starts to appear. The existence of solution where 
q = and rj ^ depends on the ratio Y/g. For instance, in Fig. 2b, one 
can find solutions for the order parameters which corresponds to the PAIR 
phase only for small values of Y. In Fig. 2c and Fig. 2d, the strength g is 
increased and solutions with rj ^ starts to appear in a larger region of the 
diagram. Therefore, the sequence Fig. 2a-Fig. 2d shows clearly that it is 
requested greater values of g to the PAIR phase occupy a larger region than 
the SG phase. The QCP given by Y = Y c is the same as g = and the SG 
phase remains unstable according to the calculated A-T line. 

Howewer, one interesting effect in the interplay between SG and the PAIR 
phase can be seen if one considers that Y and g are no longer independent 
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parameters. Assuming the following relationship: 

r = ag + T (43) 

As a justification of Eq. ()43j). one should recall the derivation of the effective 
model (see Eq given in the Appendix of the Ref. |2]. The s-ti exchange 
part, after the integration of the conducting electrons, originates the pairing 
interaction as well as the RKKY coupling between the localized spins. 

The results have been shown in Fig. (3) in a diagram T/J versus g/J. The 
position of the QCP (g c ) and the tricritical point (T tc , g tc ) is very sensitive to 
the choice of the factors a and r . These factors have been adjusted to obtain 
a second order transition between the NORMAL and PAIR phases with the 
tricritical point located in the same scale of the Figs. (1-2). Therefore, for 
a = 0.0903 and Tq = 1.9254, one can see that as long as g increases (r 
also start to increase), the results show the freezing temperature Tf being 
depressed to zero at g c . For g > g c , the solutions for the order parameters 
indicate a NORMAL phase until one gets a line transition g\(T) between 
the NORMAL and the PAIR phases and gt c > g c - This phase diagram build 
with only one independent parameter (<?) is more adequate to address the 
experiments. 

4 Conclusions 

In this work, we have investigated the competition between pairing formation 
in real space and spin glass order when tunneling is tuned by the transverse 
field T. We used the same framework of Refs. [21 El El, therefore, the parti- 
tion function is obtained using the functional integral formalism and the spins 
operators are represented by Grassmann variables. One important point is 
the use of the static approximation and the symmetry replica "ansatz" in our 
approach. It is known that treatment which neglects the dynamical behavior 
of correlation functions is not correct at low temperatures. Nevertheless, our 
interest is mainly to capture the effects that appears on the phase bound- 
aries when quantum tunneling is present due to the transverse field T. This 
procedure to find the phase boundaries is justified by the critical behaviour 
of quantum rotor model fH\ . 

The main results can be seen in Figs (1), (2) and (3). The first two figures 
show that the pairing formation is not favored when the quantum tunneling 
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is increased. At the same time, the temperature where is found the non- 
trivial spin glass ergodicity breaking decreases toward zero. For instance, in 
the Fig. 1, for the case T = where it is found a SG phase, for small g (pair- 
ing interaction strength), below the freezing temperature Tj = 0.95 J. The 
presence of transverse field (r = J) would favor the spin flipping destroying 
the double occupation of the sites. Therefore, it is necessary to increase g to 
find solutions for the order parameters where there is pairing long range or- 
der which corresponds to the PAIR phase. In that sense, the transverse field 
inhibit the pairing formation which makes the sites insensitives to a magnetic 
interaction. This results can be better seen in Fig. (2) which show clearly 
that in order to found a PAIR phase when T is increased it is necessary 
greater values of the parameter g. The position of the tricritical point found 
in the PAIR line transition is quite sensitive to the presence of the transverse 
field. It moves up with F enlarging the first order transition region in the 
phase diagram. 

In the Fig. 3, it is assumed a linear relationship between T and g (the 
strength of the pairing interaction) given in Eq (J43|) . Therefore, the strength 
of spin flipping is now related with the strength of the pairing interaction, 
he first effect when g is increased is to lead the boundary line PARA-SG to a 
QCP at g c . The pair formation is still inhibited even if g is kept increasing for 
values greater than g c . The PAIR phase only appears at the line transition 
g = gi(T). This resulting phase diagram displays phase boundaries similar 
to the experimental one for Ui_ x La x Pd,2Al 3 when x > 0.5 6j. 

In conclusion, we have studied a fermionic representation of Ising spin 
glass (SG) in the presence of transverse magnetic field T together with lo- 
cal pairing interaction. We expect that results obtained in this model can 
contribute for the study of the interplay between spin glass and supercon- 
ductivity in strongly correlated systems. Particularly, to describe the phase 
boundaries which is the main interest of this work. It should remarked that 
we have used the replica symmetry ansatz in the present work. There are 
results 0] indicating the Parisi replica permutation symmetry breaking affect 
the boundary between superconductivity and the SG phase. This is an indi- 
cation that would be necessary to go beyond the replica symmetry solution 
in the present work. That will be subject for future work. 
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